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INTRODUCTION 


The Taylor series has been generalized in many directions. It 
is our purpose here to extend’ the series to include two points about 
which a function may be expanded. The virtue of this procedure lies 
in the fact that the accuracy of using 2n derivatives in a one-point 
Taylor series is attained using only n derivatives with the two-point 
expansion. , 

The classical interpolation problem is concerned with approximating 
a fmetion^ f, with a polynomial of degree n such that the polynomial 
is in agreement with the values of the function at each of the points 
i=0;, 1^ 2, n. If in addition to the fimction values at each 

point there are available a given number of derivatives of the function 
at each pointy we have the well-known general Hermite interpolation 
formula. Thus the n+1 term Taylor series is simply a special case of 
the Hermite formula with only one point Xi and n derivatives of f 
given at xi. The special case of the Hermite formula for two points 
Xi and Xa can be written in the forin^ 


where 


and 
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This form s-'ji'ferr- from the fact that for each n the coefficients mast 


he reeaj.cu'_ated, i„e., the ■.oafficients are non-final. 

It is the purpose of this report to develop a series expansion of 
a function f about two points, xi and xb , such that the coefficients 
are final. Hence increasing the number of derivatives used involves 
only the addition of more terms to the original series. 


CALCULATION OF SKRIES COEFFICIENTS: 


d^f 


Assiime that — r , i - 0, 1, ..., k, are known at the two points 
dx^ 

Xi and X 2 . Then 


k-1 


y(x) = ^ [a^(x-X 2 ) + b^(x-xi)] [(x-xi) (x-X 2 )]^ 


( 1 ) 


i=0 


is a (2k-l.)th degree polynomial approximation of f(x) if a^ and b^ are 

d^ d^f 

determined so that — ? and — r agree at Xi and X 2 for i = 0, 1, . . . , 

dx^ dx.^ 


k-1. 


For ease in the calculations which follow, we introduce the 
transformation 


Then 


and 


^ letting X2-X1 = h. 

x-xi = ht 
X-X2 = h(t-l), 

k-1 

y(t) = y [A.(t-l) + B.t] [t(t-l)f 

1 

i=0 


( 2 ) 
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The coefficients of the transformed series (2) are related to those of 
the original series (l) by 


A. 

1 


= h'2i^l 


a, 

1 


B. = 


b. . 
1 


( 5 ) 


This transformation effectively transforms the interval fxi, X 2 ] onto 

,0, 


[ 0 , 1 ]. Now 


^ ^ dx 

dt dx dt ~ 


dx 


and in general 


o; = 

df^ 


= h‘ 


dx*^ 


Rearranging ( 2 ) in the form 


k-1 


k-1 


y(t) = ^ t^(t-l)^'^^ + ^ B^.t^'^^(t-l)^ 

i=0 i=0 

and differentiating we obtain 


k-1 


d t^ 


i=0 


Ad = ^ A, ^ + y B, id 


k-1 


i=0 


J 

dt‘ 


It can be shown from the binomial theorem that 


[t^-'^(t-l)^] 


dt" 


= 0-i-i) 


t=o 


when 2i+l ^ j ^ i+1 


^ [t^(t-l)^^^] 


dt" 




t=o 
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■when 2:.+l ^ 4 ^ 



t-i 



when 2i+l Sr j ^ i 





tr--l 



when 2i+l ^ J ^ i+1 . 


These expressions vanish when j falls outside of the given range. 
Hence we have 
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from which 
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(h) 


for j = 0, 1, 2, . . . j k 

where ^ is deleted for j=0 and [■i] denotes the integer part of ^ . 

The remainder for the classical interpolation problem involving m 

distinct data points with a.-l, j=l, 2^ m, derivatives available 

J 

2 

at each point is well-known . For our case where m = 2 and = k 

the remainder is 


R(x) = ■ [(x-xi) (x-xa)] 

( 2k) 

where 5 lies in the interval defined by x, Xi, and Xs, and f is 
continuoTis on that interval. Thus the finite form of the two-point 
Taylor series is 


k-1 


f(x) = y [a.Cx-Xa) + b.(x-xi)] [(x-xi) (x-xa)]^ 

c—j i X 

li) 


i=0 


O) 


f( 2k) I 

(2k), 


[(x-xi) (x-xa)] , 


where the a. and b. are obtained from formulas ( 5 ) s^ad (4). 
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EXAMPIiE OF APPLICATION ; 

The solution of the two -body problem, given initial coordinates and 
velocity components at a time t=to is 

X = x(t), y = y(t), z = z(t) 


5 


and for short intervals of time these solutions may be expanded in a 
Taylor series about tos 

x(t) = Xo + Xo (t-to) + ^ (t-to)^ + .. o . 

yCt) ^ Yc + Yo ft-to) + ^ (t-to)^ + ... 

z(t) zo + Zq (t-to) + ^ (■ -"'o)^ + ... 

From the differential equatio.ns representing the motion of the two 
body problem we have 


X 






3 

r- 


X = - iJLVX^ y = M*vy^ z = - jjLVz 

where v = — , and |ji is the gravitational constant times the mass 
r 

of the attracting body, 

th 

To obtain an algorithm for the n— time derivative of x, y^ and z 
we consider a power series solution of the equations of motion. 

Let 
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Substituting into the equations of motion and equating coefficients for 
the povers of (t-to) we obtain 
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and 


^ {i(X Xo + Y.Yo + Z.Zo) 


Y tx. 
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As a special case of the preceding two-body problem consider the 
motion to be confined to the x-y plane and suppose that we are given 


and 


yi^ yi Q-'t time tx 


H2.y and y2 at time t2- 


For convenience we set the period of the orbit T = 2 TT and the semi- 
major axis a = 1 so that the constant |j, = 1 . Obtaining the quantities 


d X! 


^ yi 


d XI 


dt 


n 


t=t-i 


dt 


n 


t=ti 


dt 


n 


and 


t=t2 


dt 


t=t2 


by the general procedxrre outlined above ^ we may in turn calculate the 
coefficients a^ and b^ for L>^e two-point Taylor expansion from equa- 
tions (j) and (h). 

A comparison of the resulting solutions for the two-body problem 
from the two-point Taylor series and one-point Taylor series expanded 
about ti is presented in Figure I. Here the orbit is of eccentricity 
e = ^ and ti^ t2 were chosen to be time at perigee and time at apogee 
(o and TT) respectively. The figure is a plot of the error in the cal- 
culation of the magnitude of the radius vector vs, the time of the 
orbit. All calculations were performed in double precision on an 
IBM 7094. 
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